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Performing in this paper the duality and statefinder diagnostic for a class of Quintom models 
of Dark Energy which is constructed by the spinor field, we study the possible connections among 
different Spinor Quintom models and differentiate the Quintom dark energy in spinor scenario from 
the others. A class of evolutionary trajectories of these Quintom models in Spinor field are plotted in 
the statefinder parameter planes. We also find that the power- law- like potential also plays important 
role in this class of models from the statefinder viewpoint. 
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I. INTRODUCTION 

There is mounting datum from type la supernovae 
and cosmic microwave background (CMB) radiation and 
so on [J 0, H, H[ have provided strong evidences for the 
present spatially flat and accelerated expanding universe 
which is dominated by dark sectors. Combined analysis 
of the above cosmological observations support that the 
energy of our universe is occupied by dark energy (DE) 
about 73%, dark matter about 23% and usual baryon 
matter only about 4% which can be described by the 
well known particle theory. In the context of Friedmann- 
Robertson- Walker (FRW) cosmology, this acceleration is 
attributed to an exotic form of negative pressure, the so- 
called DE. So far, the nature of dark energy remains a 
mystery. Theoretically, the obvious candidate for such a 
component is a small cosmological constant A (or vac- 
uum energy) with equation of state w = —1, while 
the cosmological model that consists of a mixture of 
Vacuum energy and cold dark matter (CDM) is called 
LCDM (or A CDM), but it is at the expense of the 
difficulties associated with the fine tuning and the co- 
incidence problems. The inspiration from inflation has 
suggested the idea that dark energy is due to a dy- 
namical component with one- or multi-scalar nelds.such 
as the Quintessence [1, 0|, PhantomQ, K-essence[l, Q. 
However, a positive kinetic term in the Lagrangian of 
Quintessence dark energy model may violate the strong 
energy condition. On the other hand, the negative ki- 
netic term in Phantom scenario leads to some quantum 
instabilities such as the violation of the dominant en- 
ergy condition, as well as the occurrence of the phe- 
nomenon of Big Rip [3. There are also other models, 
including Chaplygin g as|l l|. braneworld models [H[l3, 
holographic models (l4l. Ilot lid [l7|. and so on, which are 
resorted to accounting for the present cosmic accelerating 
expansion. 

Although the recent fits to the data in combination 
of WMAPpJ, [H, the recently released 182 SNIa Gold 
sample [ld| and also other cosmological observational 
data show remarkably the consistence of the cosmological 
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constant, it is worth of noting that a class of dynamical 
models with the equation-of-state (EoS) across —1 Quin- 
tom is mildly favored [U [22|, [H, HH . In the literature 
there have been a lot of theoretical studies of Quintom- 
like models. For example, motivated from string theory, 
the authors of Ref. |25| realized a Quintom scenario by 
considering the non-perturbative effects of a generalized 
DBI action. Moreover, a No-Go theorem has been proved 
to constrain the model building of Quintom [26j], and ac- 
cording to this No-Go theorem there are models which in- 
volve higher derivative terms for a single scalar field [27| , 
models with vector field [HI, making use of an extended 
theory of gravity f29|. n on-local string field theory [30l| . 
and others (see e.g. [3l|, M, HI H M, HI Hi], [M 
The similar work applied in scalar-tensor theory is also 
studied in Ref. 0|. 

Previously, it has been considered that a Quintom dark 
energy model and its combination with Chaplygin gas 
fluid can be realized with non-regular spinor matter [4 1|. 
Interestingly, we find that this type of model can real- 
ize many kinds of Quintom scenario with transforming 
the form of potential of the Spinor. To understand the 
possible combinations among different types of Quintom 
model in spinor field, in this paper we study the implica- 
tions of cosmic duality with this class of models. Cosmic 
duality is a mathematic feature which origins from string 
cosmologyji^, [HI and was later considered to link the 
standard cosmolo gy and phantom cosmology together, 
see Refs. [S IE S|. In Ref. [H, it was pointed 
out that there is a behavior of dual between the models 
of Quintom- A and Quintom-B. By studying the behavior 
of the energy density and pressure in spinor field we find 
a dual of the Quintom-A and Quintom-B. Meanwhile, 
we realize additional Quintom models by the aid of this 
property. Since more and more dark energy models have 
been proposed to explain the current cosmic acceleration, 
a method to discriminate between the various contenders 
in a model independent manner was introduced by Sahni 
in jH,|4^|. The new cosmological diagnostic pair {r, s}, 
called the statefinder, is a geometrical diagnostic which 
is algebraically related to the higher derivatives of the 
scale factor a with respect to time and is a natural next 
step beyond the Hubble parameter H = - and the de- 
celeration parameter q which depends on a. In this way, 
the dark energy can be distinguished more universally 
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than the model-dependent physical variables depending 
on the properties of physical fields describing the dark 
energy. We plot the trajectories in the r — s plane cor- 
responding to these kind of Quintom dark energy model. 
Departure from the fixed point {r, s} = {1,0}, for one 
given Quintom model in spinor field, in correspondence 
with the spatially flat LCDM scenario provides a nice 
way of setting up the distance from LCDM. 

This letter is organized as follows: In section 2, we in- 
vestigate the implications of the cosmic duality in Spinor 
Quintom models of dark energy. In section 3, we apply 
the statefindcr diagnostic to the Spinor Quintom dark 
energy models. Section 4 contains discussions and con- 
clusions. 



II. DUALITY OF SPINOR QUINTOM 
UNIVERSES 



To begin with the discussion, we consider a uni- 
verse filled with Quintom dark energy in spinor field (Hoi. 
[HTI . I52I ] , neglecting the contributions of the components 
of matter and radiation. And we deal with the ho- 
mogeneous and isotropic Friedmann-Robertson- Walker 
(FRW) space-time, assuming the space-time metric as 



ds 1 



dt 2 - a 2 (t)dx 2 



(1) 



In the aid of the dynamics of a spinor field which is min- 
imally coupled to Einstein's gravity [H, [H, [55[, we can 
write down the following Dirac action in a curved space- 
time background 



S v , = d A xe [-{QVD^ - D^FV) - V] 



(2) 



Here, e is the determinant of the vierbein e° and V stands 
for any scalar function of ip, ip and possibly additional 
matter fields. We will assume that V only depends on 
the scalar bilinear ij)ip. For a gauge-transformed homoge- 
neous and a space-independent spinor field, the equation 
of motion of spinor reads 



In this section we study the duality of the Quintom uni- 
verse with spinor matter 1 . Firstly, with the Lagrangian 
given by Eq. (|2|), we review the connection between 
Quintom-A and Quintom-B with the duality which has 
been studied by Ref. [33[ with two fields consists of 
quintessence-like and phantom-like. Here, we give the 
new dualities between the two models. Furthermore, the 
duality of Quintom model will be extended to the Spinor 
Quintom twice crossing —1. 

In the framework of FRW cosmology, the Friedmann 
equation reads 



(0) 



where we use units 8ttG = h = c = 1 and all parameters 
are normalized by M p — 1/\/8ttG in the letter. 

Following the work of Ref. [5?], [58| there is a form- 
invariant transformation by defining a group of quantities 
as follow: 



P 
H 



P(p)_ > 



(7) 
(8) 



Under this transformation, we obtain the corresponding 
changes for the pressure p and the EoS w, 



w 



P dp 
P dp 



(9) 
(10) 



Taking p = p in Eqs. (J9]) and (fT0|) as an example of 
detailed discussion without loss of the generality of the 
physical conclusion and information, we can get the dual 
transformation: 



H = -H, 
P = -2p-p 

w = —2 — w = — 1 



V'^ip - V, 
V'4>ijj 
V ' 



(11) 
(12) 

(13) 



Consequently, the dual form of Lagrangian reads, 



ip + -Hi/j + ij°V'ip = 0, 
■ijj + - ij°V'ip = 0, 



(3) 
(4) 



where a dot denotes a time derivative while a prime de- 
notes a derivative with respect to tpip, and H is Hubble 
parameter. To take a further derivative, we can obtain 
the solution of equation of motion: 



ipip 



N 

^3 



(5) 



(14) 



Contrasting the Lagrangian derived from Eq. @ and 
its dual form in (|14p . we may see, with the dual trans- 
formation, if the original Lagrangian is for a Quintom-A 
model the dual one is for a Quintom-B one, and vice 
versa. With this property, it is possible that the early 
Universe can be linked to other epoches of the universe. 
While from the Eq. (JT3J) , one can expect a symmetrical 



where iV* is a positive time-independent constant and we 
define it as present value of t/^. 



1 We notice such a duality has been investigated in Ref. [56| . 
however, in phantom cosmology. 
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evolutions of the EoS comparing with its dual: 

(i) . there is 

V < -► V > , 

which gives a Quintom-A scenario by describing the uni- 
verse evolves from Quintessence-like phase with > — 1 
to Phantom- like phase with < — 1; 

(ii) . there is 

v > o -> y < o , 

which gives a Quintom-B scenario for which the EoS is 
arranged to change from below —1 to above — 1; 

For our detailed discussions, we will consider three kind 
of special forms of power-law-like potentials and perform 
its semi-analytic solution to present the dual characteris- 
tics, then we take special form of potentials to study its 
numerical solutions. 

In the first instance, we take the potential as V = 
Vo{(ipil> — b) 2 + c] which can realize Quintom-A scenario, 
and the detailed discussion can be found in Ref. [4l| . Its 
dual form of the solution is a description of universe in the 
case of Quintom-B. According to Eq. ([5]), one finds that 
Tpip is decreasing along with an increasing scale factor a 
during the expansion of the universe. From the formula 
of V, we deduce that at the beginning of the evolution 
the scale factor a is very small, and so ipip becomes very 
large and make sure V' > at the beginning. Then -ipip 
decreases along with the expanding of a. At the moment 
of xpij} — b, one can see that V = which results in 
the EoS = —1. After that V' becomes less than 
0, so the universe enters a Phantom-like phase. Finally 
the universe approaches to a de-Sitter space-time in the 
Quintessence Phase in the future. Accordingly, the EoS 
of the dual form evolves from the w < — 1 and crosses -1 
by t — > from below —1 to above —1 then increases and 
declines finally approaches to the cosmological constant 
when t — > +oo. It is shown that either Quintom-A or 
Quintom-B will avoid a big rip when w < — 1. In Fig. [T] 
we plot the concrete picture of this dual pair. One can 
find this two models dual to each other rigorously. 

In succession, if V — Vq[— (ijj'ip — btyip + c], one can 
obtain a Quintom-B model (see Ref. [4l|). Take its dual 
form theoretically given by the above transformation, we 
can present the numerical solution in Fig. O Clearly, 
the duality of this Spinor Quintom model show the evo- 
lutionary picture of Quintom-A. 

These two class of models describe different behaviors 
of the cosmological evolution with one in the expanding 
phase while the other in the contracting one depending on 
the potential and initial conditions we choose. It is found 
that Quintom model and its dual form are symmetric 
around w = — 1. 

For an extending investigation, we take V = Vo[(ipijj — 
b) 2 ^[)ip + c], which can realize a picture across w = — 1 
twice, for discussion. In Fig. [3J one can see its dual model 
that evolves from its symmetrical side with respect to -1 
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FIG. 1: Plot of the evolution of the EoS of Quintom- A and its 
dual Quintom-B as a function of time for V = Vo[{ipip — b) 2 +c]. 
In the numerical calculation we take Vo = 1.0909 x 10 -117 . 
For the model parameters we choose b = 0.05, c = 10~ 3 . For 
the initial conditions we take (ipip)o = 0.051. 
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FIG. 2: Plot of the evolution of the EoS of Quintom-B and 
its dual Quintom-A as a function of time for V = Vo[— (tpip — 
b)ifnp + c]. In the numerical calculation we take Vo = 1.0909 x 
10 -117 . For the model parameters we choose b = 0.05, c = 
10 -3 . For the initial conditions we take (4>ip)o = 0.051. 



and crosses -1 two times and ultimately approaches to 
the cosmological constant boundary when t — > +oo. 

From the above analysis, we may investigate the con- 
nections among different periods or phases of our universe 
with the help of this character. As is known that with dif- 
ferent type of dark energy, the fate of the universe will be 
different. Our study in this section helps understand the 
properties of various dark energy models and their con- 
nections to the evolution and the fate of the Universe. 
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FIG. 3: Plot of the evolution of the EoS of NEC satisfied 
Quintom model as a function of time for V = Vo[(ipip — 
b) 2 ipil) + c] . In the numerical calculation we take Vo = 1.0909 x 
10~ 117 . For the model parameters we choose b = 0.05, c = 



1CT 



For the initial conditions we take (^/")o = 0.051. 



Moreover, the past and future properties can be under- 
stood by studying the above characters. One application 
of combining these properties together is Quintom-like 
bouncing cosmology, which has been intensively studied 
in Refs. [H, H|| |37[. In the meanwhile, we also realize 
three additional Quitom models. 



III. STATEFINDER DIAGNOSTIC TO SPINOR 
QUINTOM MODELS 

Based on the above discussions, it can be seen that 
there are so many dark energy models proposed to ex- 
plain the cosmic acceleration, thus how to secern these 
models become a widely attentional issue. With this re- 
gards, Sahni proposed a geometrical-constructed from 
space-time metric directly- statefinder diagnostic pair 
{r, s} , which is defined as[H,|4§| 



aH 3: 



r-1 



(15) 



Here q is the deceleration parameter with the definition 



9 = 



aH 2 ' 



(16) 



Accordingly, by showing different evolutionary trajecto- 
ries qualitatively in the r — s and r — q diagram this 
statefinder pair can distinguish different dark energy 
models from the others. Hitherto, some dark energy 
models, such as Quintessence, Phantom, the Chaply- 
gin gas, braneworld models, holographic models, and in- 
teracting and coupling dark energy models, have been 



perfectly differentiated, correlative works have been per- 
formed by Ref . [U, HI HD, M, HI HI, HI HI HE HI HI 

FroL ITU - In what follows we will apply the statefinder 
diagnostic to three Quintom models in spinor field — 
Quintom- A, Quintom-B, and the Quintom model cross- 
ing — 1 two times. We will use the form of the statefinder 
parameter written by pressure and energy density in the 
following text, 



1 



9(p + p)p 
2pp ' 



+ p)p 



pp 



(17) 



where the energy density and pressure are given by Ref. 
0. 

Taking components of dark matter and dark energy 
into account in a spatially flat universe, we can write 
down the Friedmann equation: 



H 2 



1 



(fii/, + Pro)- 



(18) 



Here p m is the energy density of dark matter with EoS 
Pm = (im ~ l)pm- Ignoring the interaction between the 
two dark sectors, we can see the energy both of dark 
energy and of dark matter are conserved and satisfy its 
continuity equation respectively, 



p^ + ?>H(p^ +p^) = 0, 



p m + 3Hp m = 0. 
As a result, we obtain the following expressions, 

9 9 — 

r = 1 + + w^)Q,^ + -w ipipil^ , 

w' - 

s = 1+ wj) H 4>ip , 

and the deceleration parameter 

1 3 o 

1= 2 + 2 W ^ V " 



(19) 
(20) 

(21) 
(22) 

(23) 



where fi^ = S±. is the fraction of dark energy density, 



and w' is the derivative of equation of state with respect 
to ijnp. 

Hereinafter, we will study the statefinder diagnostic for 
the Spinor Quintom models with three different poten- 
tials. Firstly, we discuss the Quintom-A model with the 
form of potential V — Vo[(ipip — b) 2 + c], where Vo, b, c are 
undefined parameter. In Fig. 2] and Fig. O we show the 
time evolution of statefinder pair {r, s} and {r,q}. One 
can see in r — s diagram the trajectory of this case will 
be close to the LCDM fixed point. 

Next, the trajectories of Spinor Quintom-B model with 
potential V = Vq[— (rptp — b)-tpip + c] will be plotted. 
In numerical calculations, we take the same value with 
Quintom-A model. It can be seen that the evolution- 
ary graphics of Quintom-B is roughly opposite to that of 
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FIG. 4: The r-s diagram with the form of potential V — 
Vo [('ijjip — b) 2 + c] . Evolution trajectories of r(s) in the numer- 
ical calculation we take Vo = 1.0909 x 10~ 117 . For the model 
parameters we choose b = 0.05, c = 10~ 3 . For the initial 
conditions we take (if>ip)o — 0.051. 



FIG. 6: The r-s diagram with the form of potential V = 
Vo[— {fiip — b)^ptjj + c]. Evolution trajectories of r(s) in the 
numerical calculation we take Vo = 1.0909 x 10~ 117 . For the 
model parameters we choose b = 0.05, c = 10~ 3 . For the 
initial conditions we take {i>ip)o = 0.051. 
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FIG. 5: The r-q diagram with the form of potential V = 
Vo[{ijjip — b) 2 + c]. Evolution trajectories of r(q) in the nu- 
merical calculation we take Vo = 1.0909 x 10 -117 . For the 
model parameters we choose b = 0.05, c = 10 -3 . For the 
initial conditions we take (tpip)o = 0.051. 



FIG. 7: The r-q diagram with the form of potential V = 
Vo[— ('ipip — b)%p%p + c]. Evolution trajectories of r(q) in the 
numerical calculation we take Vo = 1.0909 x 10" 117 . For the 
model parameters we choose b = 0.05, c = 10 -3 . For the 
initial conditions we take (^V)o = 0.051. 



Quintom-A in both {r, s} and {r, q} diagrams, see Fig. [B] 
and Fig. [7] for a clear cognition. 

Finally, we turn to the case of crossing the cosmological 
boundary two times. The phase portraits of {r, s} and 
{r, q} are presented in Fig. [8] and Fig. [9] respectively, 
where the values of parameters are also the same with 
those of the case of Quintom-A. 

We can see that in r — s diagram either of the portraits 



is very close to the fixed point {r, s} = {1, 0} correspond- 
ing to the LCDM with FRW cosmological model at some 
period but not passing it. In the case of Quintom-A, it is 
found that the r — s phase portrait is always in the region 
of negative s and positive r. On the contrary, most of the 
r(s) trajectories of Quintom-B lies in the opposite loca- 
tions. While the third case gives another evolutionary 
portraits. 
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In conclusion, we have investigated the dynamics of 
Spinor Quintom dark energy models by using the new 
geometrical diagnostic method — statefinder pair {r, s}. 
The Ref. [64| has applied this method to Quintom model 
and successfully differentiate this class of models with 
other dark energy models, but it is seem not useful to 
discriminate Quintom models with different potentials. 
However, as can be seen in this section, the statefinder 
diagnostic is able to differentiate different Quintom mod- 
els with diverse kinds of power-law potential in the spinor 
scenario, as well as distinguish the Spinor Quintom mod- 
els with other dark energy models. 
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FIG. 8: The r-s diagram with the form of potential V — 
Vo[(i[)ip — b) 2 tfnl} + c]. Evolution trajectories of r(s) in the 
numerical calculation we take V = 1.0909 x 10~ 117 . For the 
model parameters we choose b = 0.05, c = 10 -3 . For the 
initial conditions we take (if>ip)a = 0.051. 
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IV. CONCLUSION AND DISCUSSIONS 

To summarize, since more and more Quintom DE mod- 
els are proposed, we established the connections among 
these models by studying the cosmic duality which con- 
nects the two totally different scenarios of universe evolu- 
tion keeps the energy density of the Universe and Einstein 
equations unchanged, but transforms the Hubble param- 
eter. Besides, applying the new geometrical diagnostic 
method — statefinder pair {r, s} to the Spinor Quintom 
model, we differentiate different Quintom models with 
different kind of power-law potentials in the spinor sce- 
nario, and distinguish the Spinor Quintom models with 
other dark energy models, as well. 
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FIG. 9: The r-q diagram with the form of potential V = 
Vo[{ipip — b) 2 Tpip + c]. Evolution trajectories of r(q) in the 
numerical calculation we take Vo = 1.0909 x 10" 117 . For the 
model parameters we choose b = 0.05, c = 10 -3 . For the 
initial conditions we take (V>V)° = 0.051. 



[1] S. Perlmutter et al, Astrophys. J. 483, 565 (1997); Adam 
G. Riess et al., Astrophys. J. 116, 1009 (1998). 

[2] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003). 

[3] A. G. Riess et al, Astrophys. J. 607, 665 (2004). 

[4] U. Seljak et al, Phys . Rev. D 71, 103515 (2005) 
[arXiv:astro-ph/0407372| . 



[5] C. Wetterich, Nucl. Phys. B 302, 668 (1988). 
[6] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 3406 
(1988). 

[7] R. R. Caldwell, Phys . Lett. B 545, 23 (2002) 

[arXiv:astro-ph /9908168] . 
[8] C. Armendariz-Picon, V. F. Mukhanov and 



7 



[9] 
[10] 

[11] 

[12] 

[13] 

[14] 
[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 

[22] 

[23] 
[24] 
[25] 

[26] 

[27] 
[28] 

[29] 



[30] 

[31] 
[32] 

[33] 
[34] 



P. J. Steinhar dt, Phys . Rev. D 63, 103510 (2001) 
[arXiv:astro-ph/0006373] . 

T. Chiba, T. Oka be and M. Yamaguchi, P hys. Rev. D 
62, 023511 (2000) [arXiv:astrc~ph/9912463] . 
R. R. Caldwell, M. Kamionkowski and N. N. Wein- 
berg, Phys. Rev. Lett. 91, 071301 (2003) 
[arXiv:astro-ph/0302506] . 

A. Y. Kamenshchik, U. Moschella and V. Pasquier, Phys. 
Lett. B 511, 265 (2001) [arXiv:gr qc/0103004] . 
G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. 
B 485, 208 (2000) [arXiv:hep-th/0005016j . 
C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. Rev. 
D 65, 044023 (2002) [arXiv:astro- ph/01050 68] . 
W. Fischler and L. Susskind, arXiv:hep-th/9806039 
A. G. Cohen, D. B. K aplan arid A. E. Nelson, PhysTRev. 
Lett. 82, 4971 (1999) [arXiv:hep-th/9803132| . 
S. D. H. Hsu, Ph ys. Lett. B 594, 13 (2004) 
arXiv:hep-th/0403052| . 
M. Li, Phys. Lett. 



B 634, 101 (2006) [arXiv:astro-p h/0407432 ; H. Wei 
and R. G. Cai, Phys. Rev. D 72, 123507 (2005) 



arXiv:astro-ph/0509328 ; X. Zhang, Phys. Rev. D 74, 



B 603, 



arXiv:hep-th/0403127 



(2004) 



D. N. Spergel et al. [WMAP Collaboration], Astrophys. 
J. Suppl. 170, 377 (2007) [arXiv: astro-ph /0603449] . 

E. Komatsu et al. [WMAP Collaboration], 
larXiv:0803.0547l [ astro-ph], 

A. G. Riess et al., arXiv:astro-ph/0611572 

B. Feng, X. L. Wang and X. M. Zhang , Phys. Lett. B 
607, 35 (2005) [arXiv:astro-ph/0404224] . 

G. B. Zhao, J. Q. Xia, H. Li, C. Tao, J. M. Virey, 
Z. H. Zhu and X . Zhang, Phys. Lett. B 648, 8 (2007) 
|arXiv:astro-ph/0612728] . 

G. B. Zhao, J. Q. Xia, B. Fe ng and X. Zhang, Int. J . 
Mod. Phys. D 16, 1229 (2007) [arXiv:astro-ph/060362T] . 
Y. Wang and P. Mukherjee, Astrophys. J. 650, 1 (2006) 



arXiv:astro -ph/0604051 . 

Z. Li, J. X. Lu, Y. 
Phys. Lett. 



S. Piao, 
B 651, 



T. 



T. Qiu 

(2007) 



G. 
17, 



B. Zhao and 
1229 (2008) 



Y. F. Cai, M. Z. Li 
and X. M. Zhang, 
|arXiv:hep-th/0701016| . 
J. Q. Xia, Y. F. Cai, T. T. Qiu, 
X. Zhang, Int. J. Mod . Phys. D 
|arXiv:astro-ph/0703202] . 
M. Z. Li, B. Feng and X. M Zhang, JCAP 0512, 002 

(2005) [arXiv:hep-ph/0503268l . 

C. Armendariz-Picon, JCAP 0407, 007 (2004) 
|arXiv:astro-ph/0405267] ; H. Wei and R. G. Cai, P hys. 
Rev. D 73, 083002 (2006) [arXiv:astro-ph/0603052] . 
R. G. Cai, H. S. Zhang and A. Wang, Commun. Theor. 
Phys. 44, 948 (2005) [arXiv:hep-th/0505186l ; P. S. Apos- 
tolopou los and N. Tetradis, P hys. Rev. D 74, 064021 

(2006) (arXiv:hep-th/0604014|; K. Bamba, C Q. Ceng, 
S. Noiiri and S. D. Odintsov. larXiv:0810.4296 1 [hep-th]. 
I. Y. Aref'eva, A. S. Koshelev and S. Y. Vernov, Phys. 
Rev. D 72, 064017 (2005) [arXiv:astro-ph/0507067] ; 
S. Y. Vernov, |arXiv:astro-ph/0612487| A. S. Koshelev , 
JHEP 0704, 029 (2007) [arXiv:hep-th/0701103 
Y. F. Cai and W. Xue, larXiv:0809.4134l [hep-th]~ 

Guo, Y. S. Piao, X. M. Zhang and Y. Z. Zhang 
Lett. B608, 177 (2005) [arXiv:astro-ph/0410654] . 
Zhang, H. Li, Y. S. Piao and X. M. Zhang, Mod . 
Lett. A 21, 231 (2006) [arXiv:astro-ph/0501652] 
Guo, Y. S. Piao, X. Zha ng and Y. Z. Zhang, P hys 
D 74, 127304 (2006) [arXiv:astro-ph/0608165 . 



F. 
Z. K 

Phys 
X. F. 
Phys 
Z. K. 
Rev. 



103505 (2006) arXiv:astro-ph/0609699 . 
[35] Y. F. Cai, T. Qiu, Y. S. Piao, M. Li and X. Zhang, JHEP 

0710, 071 (2007) [arXiv:0704. 10901 [gr-qc]] . 
[36] quintom is found to be able to give a bouncing so- 
lution and its perturbations combine aspects of singu- 
lar and nonsingular bounce models, see for example: 
Y. F. Cai, T. Qiu, R. Brandenberger, Y. S. Piao and 
X. Zhang, JCAP 0803, 013 (2008) |arXiv:0711.2"T87l 
[hep-th]]; Y. F. Cai, T. t. Qiu, R. Brandenberger and 
X. m. Zhang. larXiv:0810.4677l [hep-th] . 
[37] Y. F. Cai, T T. Qiu, J. Q. Xia and X. Zhang, 
larXiv:08 08.0819 [astro-ph]; Y. F. Cai and X. Zhang, 
arXiv:0808.2551 [astro-ph]. 
[38] W. Zh ao, and Y. Zhang, Phy s. Rev. D 73, 123509 
(2006) arXiv:astro-ph /0604460 ; H. Mohseni Sadjadi 
and M. Alimohammad i, Phys. Rev. D 74, 043506 (2006) 
[arXiv:gr-qc/0605143j ; E. O. Kah ya and V. K. Onemli , 
Phys. Rev. D 76, 043512 (2007) [arXiv:gr-qc/0612026] ; 
Y. F. Cai and Y. S. Piao, Phys. Lett. B 657, 1 (2007) 
arXiv:gr-qc/0701114 . R. Lazkoz, G. Leon and I. Quiros, 
Phys. Lett. B 649, 103 (200 7) [arXiv:astro ph/0701353] ; 
H. Zhang and Z. H. Zhu, larXiv:0704.3121l [astro-p h]; 
T. Qiu, Y. F. Cai and X. M. Zhang, [arXTv:0710.0115l [gr- 
qc]; M. R. Setare, J. Sadeghi and A. Banijamali, Phys. 
Lett. B 669, 9 (2008) [arXiv:0807.0"077l [hep-th]]. 
[39] H. H. Xiong, T. Qiu, Y. F. Cai and X. Zhang, 
larXiv:0711.4469l [hep-th]; H. H. Xiong, Y. F. Cai, T. Qiu, 
Y. S. Piao and X. Zh ang, Phys. Lett. B 666, 212 
(2008) [arXiv:0805.04l3l [astro-p h]]; S. Li, Y. F. Cai 
and Y. S. Piao, arXiv:080 6T2363l [hep-p h]; S. Zhang and 
B. Chen, Phys. Lett. B 669, 4 (2008) |arXiv:0806.4435l 
[hep-ph]]. 

E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 
70, 043539 (2004) [arXiv:hep-th/0405034| . 
Y. F. Cai and J. Wang, Class. Quant. Gray. 25, 165014 
(2008) |arXiv:0806.3890l [hep-th]]. 
G. Veneziano, Phys. Lett. B 265, 287 (1991). 
J. E. Lidsey, D. Wands and E. J. Copeland, Phys. Rept. 
337, 343 (2000) [arXiv:hep-th/9909061j . 
L. P. Chimento and R. Lazkoz, Phys. Rev. Lett. 91, 
211301 (2003) [arXiv:gr-qc/0307111| . 
M. P. Dabrowski, T. Stachowiak and M. Szydlowski, 
Phys. Rev. D 68, 103519 (2003) [arXiv:hep-th/0307128] . 
L. P. C himento and D. Pavo n, Phys. Rev. D 73, 063511 
(2006) jarXiv:gr-qc/0505096 . 



Y. F. Cai, H. Li, Y S. Piao and X. M. Zh ang, Phys. Lett. 

B 646, 141 (2007) arXiv:gr-qc/0609039| . 

B. Feng, M. Li, Y. S. Piao and X. Zhang, Phys. Lett. 



[40 

[41 

[42 
[43 

[44 

[4& 

[46 

[47 

[48 

[49 

[50 
[51 
[52 
[53 
[54 



M. P. Dabrowski, C. Kiefer and B. Sandhofer, Phys. Rev. 

D 74, 044022 (2006) [arXiv:hep-th/0605229| . 

T. Chiba a nd T. Nakamura, Prog. Theor. Phys. 100, 

1077 (1998) [arXiv:astro-ph/9808022] . 

V. Sahni, T. D. Saini, A. A. Starobinsky and U. Alam, 

JETP Lett. 77, 201 (2003 ) [Pisma Zh. Eksp. Teor. Fiz. 

77, 249 (2003)] [arXiv:astro-ph/0201498] . 

C. Armendariz-Picon and P. B. Greene, Gen. Rel. Grav. 

35, 1637 (2003) [arXiv:hep-th/0301129l . 

B. Vakili, S. Jalalzadeh and H. R . Sepangi, JCAP 0505, 

006 (2005) [arXiv:gr-qc/0502076l . 

M. O. Ribas, F. P. Devec chi and G. M. Kreme r, Phys. 
Rev. D 72, 123502 (2005) [arXiv:gr-qc/0511099l . 
S. Weinberg, Gravitation and Cosmology, Cambridge 
University Press (1972). 

N. Birrell and P. Davies, Quantum Fields in Curved 



8 



Space, Cambridge University Press (1982). 
[55] M. Green, J. Schwarz, E. Witten, Superstring Theory 

Vol. 2, Chapter 12, Cambridge University Press (1987). 
[56] L. P. Chimento, F. P. Devecchi, M. Forte and 

G. M. Kremer, C lass. Quant. Grav. 25, 085007 (2008) 

|arXiv:0707.4455l [gr-qc] ] . 
[57] L. P. Chimento, Phys. Rev. D 65, 063517 (2002). 
[58] J. M. Aguirregabiria, L. P. Chimento, A. S. Jakubi 

and R. Lazkoz, Ph ys. Rev. D 67, 083518 (2003) 

[arXiv:gr-qc/03030i0] . 
[59] W. Zi mdahl and D. Pavon Gen. Rel. Grav. 36, 1483 

(2004) [arXiv:gr-qc/0311067| . 
[60] U. Alam, V. Salmi, T. D. Saini and A. A. Starobin- 

sky, Mon. Not. Roy. A stron. Soc. 344, 1057 (2003) 

[arXiv:astro-ph/0303009| . 
[61] V. Gorini, A. Kamenshchik and U. Moschella, Phys. Rev. 

D 67, 063509 (2003) [arXiv:astro-ph/0209395] . 
[62] X. Zhang, Phys. Lett. B 611, 1 (2005) 

[arXiv:astro-ph/0503075] . 



[63] X. Zhang, Int. J. Mod. Phys. D 14, 1597 (2005) 

[arXiv:astro-ph /0504586] . 
[64] P. X. Wu and H. W. Yu, Int . J. Mod. Phys. D 14, 1873 

(2005) [arXiv:gr-qc/0509036| . 

[65] X. Zha ng, F. Q. Wu and J. Z hang, JCAP 0601, 003 

(2006) [arXiv:astro-ph/041122l] . 

[66] M. R. Setare, J. Zhang and X . Zhang, JCAP 0703, 007 

(2007) [arXiv:gr-qc/0611084] . 
[67] W. Zhao. larXiv:0711.2319l [gr-qc], 

[68] G. Panotopoulos , Nucl. Phys. B 796, 66 (2008) 

[arXiv:0712.TT77l [astro-ph] ] . 
[69] Z. G. Huang, X. M. Song, H. Q. L u and W. Fang, A stro- 

phys. Space Sci. 315, 175 (2008) |arXiv:0802.2320 i [hep- 

th]]. 

[70] Y. Shao and H. Zhong, Mod. Phys. Lett. A 23, 879 

(2008) . 

[71] C. J. Feng, larXiv:0809.2502l [hep-th]; C. J. Feng, 
larXiv:0810.2594l [hep-th]. 



